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1.  Objective 


To  aid  in  a  data-reduction  process,  an  algorithm  was  generated  to  calculate  on  a 
square  mesh  (elements  with  sides  of  length  2m)  the  area  of  overlap  for  individual 
elements  with  a  square  box  (sides  of  length  2n)  that  is  not  rotated  relative  to  the 
mesh  subject  to  the  constraint  n  <  2m. 

2.  Algorithm 

Ignoring  edge  elements,  we  can  assume  the  following  construction  (Fig.  1)  without 
loss  of  generality:  1)  a  3  x  3  mesh  with  the  coordinate  system’s  origin  at  the  center 
of  the  middle  element  and  2)  a  square  box  with  center  (x,y),  which  also  lies  in  the 
middle  element.  The  constraint  ensures  that  the  box  lies  entirely  within  the  3x3 
mesh.  We  label  the  elements  (i,j)  with  i,j  =  1,2,3.  We  define  the  overlap  area  of  an 
individual  mesh  element  (i,j)  with  the  box  as  Ay.  Since  the  box  is  not  rotated  relative 
to  the  mesh,  notice  that  the  overlap  areas  are  rectangles  with  sides  Ax;  and  Ayj  with 

Ay  =  AxAyj .  (1) 
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Fig.  1  Construction  geometry 
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The  terms  Axi  and  Ayo  will  always  be  nonzero,  but  we  need  to  check  if  Axi,  Ayi, 
AX3,  and  Ay3  are  nonzero  as  well.  The  terms  Ax;  and  Ayj  are  determined  using  Table 
1  and  Eqs.  2  and  3. 


Table  1  Construction  definitions 


If 

Then 

Else 

x  -  n  <  -m 

Axi  =  n  -  x  -  m;  x“  =  -m 

Axi  =  0;  x“  =  x  -  n 

x  +  n  >  m 

AX3  =  n  +  x-m;x+  =  m 

Ax3  =  0;  x+  =  x  +  n 

y  -  n  <  -m 

Ay3  =  n  -  y  -  m;  y“  =  -m 

Ay3  =  0;  y"  =  y  -  n 

y  +  n  >  m 

Ayi  =  n  +  y-m;y+  =  m 

Ayi  =  0;  y+  =  y  +  n 

Ax2  =  X+  -  X  . 

(2) 
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